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ABSTRACT 



o . 

i Context. Theoretical stellar evolutionary models are still affected by not negligible uncertainties due to the errors in the adopted 

physical inputs. 

Aims. In this paper, using our updated stellar evolutionary code, we quantitatively evaluate the effects of the uncertainties in the 
main physical inputs on the evolutionary characteristics of low mass stars, and thus of old stellar clusters, from the main sequence to 
the zero age horizontal branch (ZAHB). To this aim we calculated more than 3000 stellar tracks and isochrones, with updated solar 
mixture, by changing the following physical inputs within their current range of uncertainty: 'H(p,ve + ) 2 H, 14 N(p,y) 15 0, and triple-ff 
reaction rates, radiative and conductive opacities, neutrino energy losses, and microscopic diffusion velocities. 
Methods. The analysis was conducted performing a systematic variation on a fixed grid, in a way to obtain a full crossing of the 
perturbed input values. The effect of the variations of the chosen physical inputs on relevant stellar evolutionary features, such as the 
turn-off luminosity, the central hydrogen exhaustion time, the red-giant branch tip luminosity, the helium core mass, and the ZAHB 
luminosity in the RR Lyrae region are analyzed in a statistical way. 



Results. We find that, for a 0.9 Mq model, the cumulative uncertainty on the turn-off, the red-giant branch tip, and the ZAHB 
luminosities accounts for ±0.02 dex, ±0.03 dex, and ±0.045 dex respectively, while the central hydrogen exhaustion time varies of 
about ±0.7 Gyr. For all examined features the most relevant effect is due to the radiative opacities uncertainty; for the later evolutionary 
stages the second most important effect is due to the triple-cf reaction rate uncertainty. For an isochrone of 12 Gyr, we find that the 
isochrone turn-off log luminosity varies of ±0.013 dex, the mass at the isochrone turn-off varies of ±0.015 Mq, and the difference 
between ZAHB and turn-off log-luminosity varies of ±0.05 dex. The effect of the physical uncertainty affecting the age inferred from 
, turn-off luminosity and from the vertical method are of ± 0.375 Gyr and ± 1.25 Gyr respectively. 

Key words, methods: statistical - stars: evolution - stars: horizontal-branch - stars: interiors - stars: low-mass - stars: Hertzsprung- 
Russell and C-M diagrams 
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' 1. Introduction by different authors adopting different input physics and/or pre- 

U^. , scriptions for the treatment of processes occurring in stars. These 

^ . An evaluation of the global uncertainty in stellar models due to models m in fact the result of comp i ex calculations relying on 

■ the micro-physical inputs is essential for understanding the ac- many physical assumptions (i.e. equation of state, radiative and 

; tual significance of the application of these models when denv- conductive opacities, nuclear reaction cross sections, neutrino 

, ing quantitatively fundamental stellar, and even cosmological, emission rates, etc.), algorithms describing physical processes 

7^ ' parameters. In fact, theoretical models for the structure and evo- (i e convec tive transport, rotation, etc.), and input parameters 

;> . lution of stars are indispensable tools in many astronomical re- (i e inMal me tallicity, initial helium abundance, heavy-element 

search areas. Much fundamental information on resolved stellar mixture , e tc), each affected by uncertainties, 

populations otherwise inaccessible, as for example the age, is A first estimate of the unC ertainty in stellar models can be 

H , obtained by comparing observational data and theoretical pre- obtained by compar i ng the results pr ovided b y differen t au- 

dictions. Furthermore, evolutionary models play a crucial role thors and codefj (see e g Berteln et al 2009; Dotter et .,, 2007; 

also in the studies of unresolved stellar populations, since they |p ietrinferni et a fl2006h . However, this kind of approach does not 

are a fundamental ingredient for the required stellar population allow tQ disentang i e and quanti f y the contributions of the various 

synthesis tools. As such, in the last decades a huge effort has uncertainty sources 

been focused on refining the accuracy and reliability of evolu- A more suitable approach consists in changing a give n 
tionary predictions. As a result, stellar evolution theory has be- ■ h ics at a time ke j ^ the Q±eT j and 
come one of the most robust area of astrophysical research al- parameters fixed _ For hydrogen-burning models of metal 
lowing a firm understanding of the main stellar population char- poor low . mass sta rs, an early example of this technique 
actenstics. is provided by Ichabover etalJ (1 1995b and subsequently ex- 
However, the current generation of stellar models are still tended to more advanced evolution ar y phases by variou s 
affected by not negligible uncertainties, as proved by the dis- authors such as ICassisi et alJ (fl998h : iBrocato et all C998); 
crepancies among the stellar tracks and isochrones computed ICastellani & Degrinnocentill l999i. and for w hite dwarf cool- 
ing models by Prada Moroni & Straniero ( 2002). Many other pa- 
Send offprint requests to: G. Valle, valle@df.unipi.it pers followed this approach focusing on different mass ranges 
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and/or evolutionary 



phases (see e.g. Castellani 
Salaris et al.1 l2002t llmbriani 



et alJ 
etal.1 



2000; 



2004 



Imbriani et al. 2001 

Weiss et al.1 120051: iPrada Moroni & Stranierol 120071: IValle et al.l 
2009; lTognelli et al.ll201 lh .~ 

However, the previous method does not allow to quan- 
tify the possible interactions among the different input physics. 
A more systematic and exhaustive analysis consists in vary- 
ing simultaneously all the main input p hysics adopting ei 



flier a Monte Carlo technique (se e e.g. | Chabover et al.l 



19981: IChabover & Kraussl 120021: iKrauss & Cha bover 



1996, 



2003 



Biork & Chabover! or a systematic variation of the inputs 



in a fixed grid. These kind of studies are clearly much more time 
consuming than the previous ones, since they require the com- 
putation of a huge number of stellar models. This is the reason 
why a thorough analysis of this kind is still lacking. 

In the present paper we begin to fill this gap focusing only 
on the impact on stellar models of low-mass stars, from the main 
sequence (MS) to the zero age horizontal branch (ZAHB), of the 
uncertainties affecting some of main physical inputs adopted in 
modern stellar evolution codes. 

The aim is twofold: first, to give an estimate of the cumu- 
lative effect of physical uncertainties on the tracks, isochrones 
and on the main stellar quantities; second, to break down such 
a global variability, evidencing the effects of the single physical 
inputs on different stages of stellar evolution. 

Relying on the computation of a large number of stellar mod- 
els (i.e. 3159 stellar tracks and 567 isochrones), we performed 
a rigorous statistical analysis of the effects of the variations of 
the chosen physical inputs on relevant stellar evolutionary fea- 
tures, such as the turn-off luminosity, the central hydrogen ex- 
haustion time, the luminosity and helium-core mass at the Red 
Giant Branch (RGB) tip, and the ZAHB luminosity in the RR 
Lyrae region at log T e ff = 3.83. 

Section|2]is devoted to a description of the method employed 
for the calculations and of the physical inputs relevant to the cal- 
culations and of their current uncertainty. Section[3]and[4]present 
and analyze the results on the global uncertainty of our reference 
model due to all the variations of the chosen inputs within their 
estimated uncertainty. In Sec. [5] and [6] we report the correspond- 
ing uncertainty on the isochrones. Some concluding remarks are 
given in Sect. [7] 



2. Description of the method 

For the present analysis we are only interested in quantifying 
the theoretical uncertainty affecting low-mass stellar models due 
to the cumulative effects of the uncertainties in the main input 
physics. We focus on the typical member of an old stellar popu- 
lation, choosing for reference a M — 0.90 Mq model with initial 
metallicity Z = 0.006 - suitable for metal-rich galactic globu- 
lar clusters and i n the Magellan i c Clou ds - with heavy elements 
solar mixture by Asrflun d et al.l (120091) . The initial helium abun- 
dance, Y = 0.26, was obtained following the often adopted lin- 
ear helium-to-metal enrichment law given by: Y = Y p + ||Z, 
with cosmological 4 He abundance Y„ = 0. 2485 dCvburt et alJ 
12004 ISteigmanll2006l iPeimbert et lfll2007allbl) . In this work we 
assume AT/AZ = 2, a typical value for this qu antity, still affected 
by several important sources of uncertaint y ( Pagel & Portinari 
ll998H.Timenez et aUl2003blFlvndl2004tlGennaro et al.ll2010h . 

The adopted stellar evolutionary code, FRANEC, has been 
extensively described in p revious papers (ICariulo et al.l [2004: 
Degl' Innocenti et al.l [2008L and references therein). A detailed 
discussion of the recent updates of the physical inputs can be 



found in I Valle et all (120091) and lTognelli et all (l201ll) . The code 
adopted here is the same used for the construction of the Pisa 
St ellar Evolution Data BasJ] f or low-mass stars, as illustrated 
in lDelFOmodarme et alJ d2012l) . where a detailed description of 
the inputs of the stellar evolutionary code and of the ZAHB0 
construction technique can be found. 

The mixing length formalism (Bohm-Vit enselll958l) is used 
to treat the convective transport in superadiabatic regions. In 
such a scheme the convection efficiency depends on the mix- 
ing length parameter I = a m \H p , where H p is the local pressure 
scale height and a m \ a free parameter to be calibrated. 

The treatment of the superadiabatic convective transport is 
one of the weakest aspect in current generation of stellar evo- 
lution codes and consequently one of the main sources of theo- 
retical uncertainty, which, however, mainly affects effective tem- 
perature predictions while luminosity evaluations are affected in 
very much less extent. Such an uncertainty should be considered 
as systematic and due to an oversimplified treatment of a very 
complex phenomenon, whose precise physical description is not 
yet available. Since we are interested only in evaluating the cu- 
mulative effect of the main input physics uncertainties, we did 
not take into account the effect of the large error in the efficiency 
of the superadiabatic convective transport, i.e. on the a m \ value. 
For this reason, all the models computed in the present analysis 
adopt ff m i = 1.90. However, a change of the a m \ value chosen 
in the computations might in principle affect the estimate of the 
extent of the cumulative uncertainties in stellar predictions. For 
this reason, as shown in detail in Appendix ICl we computed ad- 
ditional sets of models with two different mixing-length parame- 
ter values, namely a m \ = 1.70 and 1.80, checking the robustness 
of the results presented in the paper for a reasonable change in 
the mixing-length value. 

2.1. Input physics uncertainties 

The computation of stellar models relies on the detailed knowl- 
edge of matter and radiation behavior in the temperature and 
density regimes typical of stellar interiors. These input physics 
include: radiative and conductive opacities, equation-of-state 
(EOS), nuclear reaction cross sections, neutrino emission rates 
(i.e. plasma, photo, pair, and bremsstrahlung processes). 

In this subsection, we briefly summarize the uncertainty on 
the up to date input physics, focusing only on those which play 
a relevant role in the structure and evolution of the reference 
case, with the notable exception of the EOS. In fact, in spite of 
its importance in modeling stellar interiors, a detailed and rigor- 
ous evaluation of the propagation of the EOS uncertainties into 
the final model is a difficult task. The reason is that the main 
thermodynamic quantities required for computing stellar models 
(i.e. pressure, temperature, density, specific heat, adiabatic gradi- 
ent, etc.) are available only in the form of numeric tables, where 
the values of the various quantities are given without the associ- 
ated uncertainty. Furthermore, those thermodynamic quantities 
are related in a not trivial way. Thus, we preferred to fix the EOS 
as it was without errors rather than adopting a too crude treat- 
ment. 



1 http : //astro . df . unipi . it/ stellar-models/ 

2 We followed a synthetic method to build the ZAHB models as we 
did not evolve the tracks through the He-flash. The ZAHB models have 
been computed by accreting envelopes of different mass extensions onto 
the He-core left at the tip of the RGB. 
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Present calculations used the most re cent version of the 
OPAL EOS3, namely th e 2006 release dRogers et al.1 Il996t 
iRogers & Navfonovll2002b . 

Anyway an estimate of the impact of two different choices 
of EOS is presented in Appendix iDl where we evaluate the dif- 
ferences i n the evolu tionary features computed with OPAL and 
FreeEOS (Ilrwinll2004 . 

Nuclear reaction rates Regarding the nuclear reactions rele- 
vant for hydrogen and helium burning phases, we adopted the 
same cross sections detailed in DeH'Omodar me et a DdaOll. In 
order to avoid the huge and useless explosion of computational 
time that would result taking into account the uncertainty in the 
cross sections of all the nuclear reactions belonging to the hydro- 
gen and helium burning networks, we limited the analysis only to 
the three main reactions, namely the 'H(p,ve + ) 2 H, 14 N(p, y) 15 0, 
and triple-a. In fact the 'H(p,ve + ) 2 H and H N(p, y) 15 6 are the 
lowest cross sections which thus drive, respectively, the effi- 
ciency of the proton-proton chain and the CNO cycle, while the 
triple-a- cross section influences both the core helium flash and 
the ZAHB luminosity. 

Experimental measurements of the 'H(p,ve + ) 2 H cross sec- 
tion are not available, thus one has to rely only on theoretical 
models. The uncertainty in the S p n(0) factor is of the order of 1% 
dKamionkowski & Bahcalfll994tlAdelberger et al.ll 19981 1201 lb . 
However, the current version of the FRANEC, as many other 
evolutionary codes, calculates the 'H(p,ve + ) 2 H reaction rate fol- 
lowing the a nalytical approxima tion provided by the NACRE 
compilation ( Angulo et al1 ll999h . whose accuracy with respect 
to the tabulated rates is better than 3%. This is the value of un- 
certainty we adopted as it is the dominant source of error in the 
rate calculations. 

For 14 N(p, j) l5 Q reaction rate we used the recent estimate by 
llmbriani et alj ( 2005b . The uncertainty on the analytical fit given 
in that paper is about 8-h10% in the range of temperature [10 5 - 
10 8 ] K. In the ca l culatio ns we assumed an uncertainty of 10%. 

Fvnb o et all J2005) reported new measurements for the 
triple-Q' rate: however, in the temperature range of interest, 
the differ ences with respect to the most widely adopted rate 
(NACRE, Angulo et al.l ll999) are within the uncertainty eval- 
uated in the NACRE compilation, reaching a maximum of » 
20% at temperatures of about 10 s K (see e.g. lWeiss et al] |2005). 
For this reason we decided to take the error quoted by NACRE 
(20%), in the temperature range of interest, as an estimate of the 
uncertainty on the triple-a reaction rate. 

Radiative opacity The radiative opacity is one the most im- 
porta nt input for stellar model calculations. As described 
in iDelF Omodarme et al] d2012b . we implemented in our 
code the opacity table s prov ided in 2005 by OPAL0 (see 
e.g. Ilglesias & Rogers! 11996b for logT(K) > 4.5 and by 
Ferg uson et al.1 (l2005h for lower temperatures. Both opacity ta- 
bles ha ve been compu t ed ado pting the solar heavy-element mix- 
ture bv lAsplund et al.1 d2009t) . 

In spite of the crucial importance of radiative opacity in the 
computation of stellar models, the quoted tables do not contain 
the uncertainty associated to the sing le Rosseland mean opac- 
ity coefficients. Fu r therm ore, neither [iglesias & Rogers! (1 1996b 
nor lFerguson et al] (l2005h provide even a rough estimate of the 



accuracy of the results. The same occurs for the radiative opac- 
ity tables made available by the Opacity Pr oject (hereafter OP, 
ISeaton & Badnellll2004! IBadnell et aT] l2005h . Given such a sit- 
uation, a first idea of the accuracy level of the current radiative 
opacity generations can be achieved by compa ring the r esults 
provided by different and independent groups. lRosel fcOO 1 ) per- 
formed this kind of analysis for a plasma mixture very similar 
in composition, temperature, and density to that in the Sun's 
center finding an agreement between Rosseland mean opacity 
coefficients provided by different codes at the level of 5%. The 
agreement gets worse at different plasma conditions, as shown 
by comparing the results provided by OPAL and OP (here- 
after £ 1 OPAL and k® p ) over their whole validity domain, which 
are in agreement w ithin about 10-15% (Seaton & Badnell 2004; 
IBadnell et alEool . 

To further extend the investigation, we evaluated the differ- 
ences in the (log T, log R) plane between OPAL and OP radia- 
tive opacities computed for Asplu nd et alj (|2009) solar elements 
mixture. As in OPAL and OP works, we used here the variable 
R = pjT\, where p is the mass density and = 10~ 6 x T with 
T in K. FigureQ]shows the values of (k° PAL - k° p )/k° PAL in four 
contour plots, for different hydrogen abundances (i.e. X - 0.8, 
0.5, 0.2, and 0.0) and metallicity Z = 0.004. Notice that we did 
not plot the comparison for the same metallicity of our reference 
track, i.e. Z = 0.006, because the corresponding tables are not 
available neither in OPAL nor OP calculations, thus we preferred 
to use the nearest value of the metallicity for which both groups 
provide the opacity tables, i.e. Z = 0.004, rather than interpo- 
late. The computed evolution of the standard stellar model with 
M = 0.90 Mq, Z = 0.006, Y = 0.26 is superimposed to each 
panel. To give an indication of the region spanned in the plane 
by the whole stellar structure, we showed the path of four differ- 
ent mass fractions of the structure, from the center to the outer 
layers, corresponding to the 0.99974 mass fraction (respectively, 
labels from 0.00 to 1 .00 in Fig. Q). 

Looking at the panels of Fig. Q] we see that, for the selected 
model, a 5% uncertainty on the values of radiative opacities is 
adequate, so we adopted this value in the calculations. We as- 
sumed the same uncertainty also for the low-temperature r adia- 
tive opacity coefficients computed by Fergu son et alJ (2005). 



Conductive opacity In regions of stellar interiors characterized 
by electron degeneracy, as in the helium core of low-mass red 
giant stars, thermal electron conduction, elsewhere inefficient, 
contributes significantly to the energy transport. Thus, in these 
regimes it becomes an important input which affects the stellar 
model. 

For t he electron conducti on opacities k c , we adopted the re- 
sults by ICassisi et al. (l2007h . which improved the conductive 
opacities by iPotekhiri (Q999). In the mildly degenerate, weakly 
or moderately coupled plasmas, typical of the red giant helium 
cores, the dynamical plasma screen i ng ma y be non-negligible. 
This effect was studied by Lampe dl968b . but it was not in- 
cluded in either iHubbard & Lampd 11969) or in Cassisi ~eTaT] 



3 Tables available at http : //rdc . llnl . gov/EOS_20O5/ 

4 The OPAL radiative opacity can be found at the URL: 



http : / /opalopacity . llnl . gov/new . html 



(2007) calculations. Using an estimate for the correction due 
to the dynamical plasma scr eening, Potekhin su ggested (2011, 
private communication) that Cassisi et al. (2007) opacities are 
uncertain by 5% at the center and few ten percent in the outer 
regions of the helium core. However, in our case the uncertainty 
is lower than the 10% in about the 80% of the He-core mass, 
which roughly corresponds to the zone dominated by electron 
conductivity. So we adopted a 5% uncertainty as sensible choice 
for conductive opacity. 
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Neutrino emission rates At high temperature and density, sev- 
eral processes of neutrino emission are efficient (i.e. plasma, 
photo, pair, and bremsstrahlung processes), which provide addi- 
tional and relevant energy loss channels in stellar interiors, with 
the consequent significant effect on stellar structures. During the 
red giant phase for the chosen reference case, in the dense and 
not extremely hot helium core, neutrino losses are dominated by 
the plasma neutrino emission process, whose efficiency affects 
important quantities of red giant evolution, such as the helium 
core mass and the stellar luminosity at the helium-flash onset. 

We followed the fitting formulae given in (lHaftet al.lll994h 
to comput e the plasma neutrino emission r ate. As reported in lit- 
erature dltoh et alj|1996t lHaft et alj|1994l) . the accuracy of that 
approximation is better than 4% for almost every value of tem- 
perature and density and better than 5% everywhere. In the cal- 
culations, we adopted an uncertainty of 4%. 



Helium and heavy elements diffusion velocities The com- 
putation of the diffusion velocities of helium and heavy ele- 
ments in our code is performed by the routin e developed by 
Thoul et al. (1994 ]) (for more details see e.g. ICastellani et alj 



the inputs of the calculations, since it does not assume physical 
independence of the individual processes. 

An alternative approach to the proble m would require the 
use o f Monte Carlo simula t ions (s e e e.g. [ Chabover et a J 1996, 



1998; Chabo ver & Kraussl l2002t iKrauss & Chaboverl 120031 
Biork & Chabover 12006ft : this technique allows the construction 
of probabilistic confidence intervals for the most interesting evo- 
lutionary features. 

The grid technique employed is distribution-free, i.e. it needs 
only the specification of a sensible range of variation for each 
physical input, and it does not rely on explicit specification of the 
parent distributions of the physical parameters that we varied. 
The simpler grid variation technique has the additional advan- 
tage to be less computationally expensive than a Monte Carlo. 
In fact a key value to evaluate for a Monte Carlo simulation is 
the number of runs N to perform in order to have an acceptable 
coverage of the input parameter hyperspace. As an example, di- 
viding the range of variation of all the parameters in 4 zones ac- 
cording to quartiles, the hyperspace of parameters will consist on 
n — 4 7 = 16384 hypercells. In the hypothesis that all the param- 
eters have uniform distribution, the probability of having k hit 
of a cell during the Monte Carlo simulations is given by Poisson 



19971: iDeirOmoda rme et al.l l2012h . These diffusion velocities density function with mean N/n: P{k,N/n) = k N/n /(N/n)l 



are generally th ought to be accurate at 10-^15% (see e.g. 
lThouletal.ll 1994 . thus, in the present computations we adopted 
a conservative uncertainty of 15%. 

The selected physical inputs and their assumed uncertainty 
are listed in Table Q] 

2.2. Method to evaluate the cumulative uncertainties 

The complexity of stellar evolution calculations hampers an an- 
alytical evaluation of the impact of the variation of the chosen 
inputs on stellar models calculation, so that the problem must be 
addressed by direct computation of perturbed stellar models, i.e. 
models adopting physical inputs perturbed within the uncertain- 
ties. 

Several choices are available for the design of the sample to 
investigate. As already discussed in the introduction, the widely 
followed approach is to allow for the variation of one parameter 
at a time, making possible the quantification of the separate ef- 
fect of the studied inputs. This approach is however vulnerable 
in presence of interactions among the inputs, since it does not 
allow to detect a possible synergyc effect due to the variation of 
several inputs at a time. 

To avoid this weakness, we employ a more robust, even if 
much more computationally expensive, technique namely a sys- 
tematic variation of the inputs on a fixed grid. For each physical 
input, we introduced a three-values multiplier p, with value 1 .00 
for the reference case and values 1.00 + Ap/ for perturbed cases 
(A/7, is the uncertainty listed in TableQ]), which defines the range 
of variations. For each stellar track calculation, a set of multi- 
plier values (i.e. p\, . . . , p-j for the seven input physics allowed 
to vary) is chosen and kept constant during the evolution of the 
stellar structure. In order to cover the whole parameters space, 
calculations of stellar tracks were performed for a full crossing, 
i.e. each parameter value p, was crossed with all the values of 
the other parameters pj, with j + i. In this way, we computed 
stellar models for all the possible sets of multiplier values. A to- 
tal of 3 7 = 2187 tracks were then computed, with same mass, 
chemical composition and atmi. 

The technique allows the exploration of the edge of the vari- 
ability region, and it is robust in presence of interaction among 



If we require a probability a that a cell has zero hit - i.e. it is 
unexplored during the simulations - we have: P(Q,N/n) = a. 
Solving for a - 0.15, an acceptable compromise between cov- 
erage of the hyperspace and accuracy of the output, it results 
N ~ 31000 run, which is an order of magnitude bigger than our 
sample size. 



Table 1. Physical inputs perturbed in the calculations and their 
assumed uncertainty. In parentheses are defined the abbrevia- 
tions used in the following Tables. 



description 


parameter 


uncertainty 


'H(p,ye + ) 2 H reaction rate (pp) 


Pi 


3% 


14 N(p,y) l5 reaction rate ( 14 N) 


P2 


10% 


radiative opacity (k T ) 


P3 


5% 


microscopic diffusion velocities (v d ) 


P4 


15% 


triple-a reaction rate (3a) 


PS 


20% 


neutrino emission rate (y) 


Pb 


4% 


conductive opacity (k c ) 


Pi 


5% 



3. Global physical uncertainty in stellar models 

The availability of a large set of stellar models covering all 
the possible combinations of simultaneously perturbed input 
physics allow us to quantify the cumulative physical uncertainty 
in stellar models. 

The combined effect in the theoretical plane (log T^, 
log L/Lq) of the variation of all the seven physical inputs (i.e. 
pi,...,pq in Table [TJ selected for the calculations is displayed 
in Fig. |2] where the track computed with the standard inputs is 
enveloped by an error stripe constructed by considering the re- 
gion of the plane spanned by the perturbed stellar models. To the 
best of our knowledge, this is the first time that an error stripe is 
computed and plotted for theoretical stellar tracks. A detailed de- 
scription of the technique employed for the construction of the 
stripe is given in AppendixlAl 

The considerable narrowing of the error stripe in the RGB 
is due to the fact that the perturbed stellar models are disposed 
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Fig. 1. Contour plots of the difference between OPAL and OP calculations for different hydrogen abundance X and metallicity 
Z = 0.004. The temperature T is in K, while R is in g cm -3 K -3 . The coloured scale marks the values of the relative difference 
(£, 0PAL - k® F )/k® FAL . Solid line: evolutionary path of the stellar center; dashed line: evolutionary path of the 0.50 mass fraction of 
the structure; dotted line: path of the 0.95 mass fraction of the structure; dot-dashed line: path of the 0.99974 mass fraction of the 
structure, labeled as 1 .00. 



along the tracks itself and do not imply the vanishing of the un- 
certainty. This effect is best evidenced in left panel of Fig. [3] 
which displays the range in A log L/Lq, computed with respect 
to the standard track, from ZAMS to helium flash. To perform 
the comparison, the raw tracks were reduced to a set of tracks 
with the same number of homologous points. Details about the 
reduction procedure are reported in Appendix [A] We note that 
the error stripe has a nearly constant width of about 0.05 dex 
until the central hydrogen exhaustion, while the error grows to 
about 0.07 dex in the final part of RGB, from the RGB bump to 
helium flash. 

In order to compare this uncertainty with the one due to a 
variation in the initial chemical composition, Fig. [3] also shows 
stellar models computed with standard physical input but differ- 
ent Z and Y. We trace a variation in [Fe/H] of +0.07 - a typical 
value of uncertainty for a cluster of metallicity similar to our ref- 
erence case - with unchanged AT/AZ = 2, resulting in two sets: 
the first one with Z = 0.005, Y = 0.258 and the second one with 
Z = 0.007, Y = 0.262. As one can see in the Figure, the impact of 
the quoted chemical variation on the predicted stellar luminos- 
ity is essentially the same of the input physics uncertainties until 
the central hydrogen exhaustion. At later evolutionary stages the 



former effect prevails, becoming the dominant by approximately 
20% after RGB bump. 

As for evolutionary time, in the right panel of Fig.[3]we show 
the evolution of the central hydrogen abundance X c as a function 
of time (Gyr). In this case the effect of the variation of the phys- 
ical inputs is larger than the one due to the quoted changes in 
chemical abundances. For instance, the range for the X c exhaus- 
tion time due to physics variation is [9.83 - 1 1.26] Gyr, while the 
one due to chemical variation is [10.03 - 1 1.00] Gyr. 

The radius is another quantity worth to be discussed, whose 
accurate determination is important also in determining the prop- 
erties of any orbiting exoplanet. Regarding the impact of the per- 
turbed physical inputs on the predicted stellar radius, in Fig.|4]we 
show the range in AR/R, computed with respect to the standard 
track, from ZAMS to helium flash. The error stripe has an in- 
creasing width from about 0.02 at the ZAMS to about 0.04 at 
the central hydrogen exhaustion, while the width grows to about 
0.09 in the final part of RGB . 

We quantified also the global physical uncertainty in the the- 
oretical predictions of various stellar quantities of common in- 
terest, namely the turn-off luminosity, the central hydrogen ex- 
haustion time tu, the luminosity L t ; p and the helium core mass 
M" e at the RGB tip, and the ZAHB luminosity in the RR Lyrae 
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Fig. 2. Left panel: Hertzsprung-Russell (HR) diagram showing the error stripe due to the variation of all the seven analyzed physical 
inputs (i.e. p\, . . . ,p-j in Table[T} on the stellar track with M = 0.9 Mq, Z = 0.006, Y = 0.26 from pre-main sequence to helium 
flash. The narrowing of the error stripe in the RGB is due to the fact that the perturbed stellar models are disposed along the tracks 
itself. See text for details. Right panel: as in the left panel, but for the ZAHB. 
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Fig. 3. Left panel: A log L/Lq with respect to the track computed with standard inputs from ZAMS to helium flash. The comparison 
was performed on a set of reduced tracks with the same number of homologous points. Right panel: central hydrogen abundance X c 
vs. time (Gyr). In both panels, the red solid and green dashed lines show for comparison the output of models with standard physical 
inputs and Z = 0.005, 0.007 respectively. 
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Fig. 4. AR/R with respect to the track computed with standard 
inputs from ZAMS to helium flash. 



region Lhb- In the case of turn-off luminosity we chose to inves- 
tigate the luminosity of a point brighter and 100 K lower than the 
turn-off (her eafter BTO), adopting a technique similar to the one 
proposed by Chab over et al.l d 19961) for isochrones in the (B-V, 
V) plane. This method has the advantage of reducing the intrin- 
sic luminosity variation of the canonical turn-off, whose position 
is difficult to accurately identify both in observed Hertzsprung- 
Russell (HR) diagrams and in theoretical tracks/isochrones, as 
a consequence of the almost vertical slope (i.e. large luminosity 
variation at essentially the same effective temperature). Then, 
a small fluctuation in the effective temperature determination 
could have a large effect in the determination of the turn-off lu- 
minosity. 

Table |2] lists the total range of variation in predictions of the 
selected quantities for our reference stellar track due to current 
input physics uncertainties. The turn-off log luminosity log Lbto 
varies in the range [0.334 - 0.376] dex (range half- width 0.021 
dex, * 6% of the value obtained with unperturbed physical in- 
puts). The total range of variation of the predicted central hydro- 
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Table 2. Total range of variation and range half- width of the the- 
oretical predictions for the selected quantities for our reference 
case, i.e. M = 0.90 M© with Z = 0.006 and Y = 0.26, due to 
input physics uncertainties. 



quantity 


variation range 


range half-width 


log L BTO 


[0.334 - 0.376] dex 


0.021 dex 


to 


[9.83 - 11.26] Gyr 


0.72 Gyr 


log Lap 


[3.38 - 3.44] dex 


0.03 dex 


M c He 


[0.4796 - 0.4879] M Q 


0.0042 Mq 


log L HB 


[1.52 - 1.61] dex 


0.045 dex 



gen exhaustion time f H is [9.83 - 1 1.26] Gyr (0.72 Gyr, * 6.5%). 
The RGB tip log L t j p and ZAHB log Lhb log luminosities vary, 
respectively, in the ranges [3.38 - 3.44] dex (0.03 dex, « 1%) 
and [1.52-1.61] dex (0.045 dex, * 3%). Finally, the helium core 
mass at the RGB tip M^ e varies in the range [0.4796 - 0.4879] 
M Q (0.0042 Mq, « 0.85%). 

4. Statistical analysis of physical uncertainty in 
stellar models 

The large set of computed stellar models is suitable for an accu- 
rate statistical analysis of the effect of the variation of the chosen 
physical inputs on relevant stellar evolutionary features. Beside 
the quantification of the whole range of uncertainty performed 
in the previous section, it is possible to disentangle the effects of 
the different physical inputs on the above selected stellar quanti- 
ties. 

The dependence of the aforementioned evolutionary quanti- 
ties on physical inputs was explored by means of linear regres- 
sion models. These models were constructed extracting the val- 
ues of the chosen dependent variable in study (i.e. Lbto, f H, Ltip, 
M He , and Lhb) from the computed stellar tracks and regressing 
it against the independent variables (more properly defined as 
predictor variables or covariates), in our case the values of the 
parameter p t . 

The regression model construction started from models lin- 
ear in the physical inputs, since a priori we do not anticipate 
the need of higher power of the inputs. Similarly, at first no in- 
teraction among the covariates were considered. These choices 
were supported by the detailed a posteriori analyzes of each 
regression model: it results in fact that the insertion of higher 
power of the predictors or explicit interaction among them is 
not needed, because the aforementioned linear models capture 
almost the whole variation of the data. These models included 
as covariates only the physical inputs which can have an influ- 
ence on the studied evolutionary feature: for Lbto and fa only 
the first four parameters of Table [TJ while for later evolutionary 
stages all the parameters were used. The models were fitted to 
the data with a least-squares metho d using the software R 2.14.1 
dR Development Core TearrfeOl ll) . Results of the statistical an- 
alyzes were considered statistical significant for p-value < 0.05. 

For the BTO log-luminosity log Lbto and tu the regression 
models were: 



log Lbto, t H = /3 + ^yS, p t 



(1) 



i=l 



where /?o> ■ • ■ ,/?4 were the regression coefficients to be estimated 
by the fit and p\, . . . , p4 respectively the perturbation multipliers 
(see Table [TJ for 'H(p,ve + ) 2 H and 14 N(p,y) I5 reaction rates, 
radiative opacity k r , and microscopic diffusion velocities. 
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Fig. 5. Boxplot of the impact of the variation of the most im- 
portant physical inputs on log Lbto and on tu- The black thick 
lines show the median of the data set, while the box marks the 
interquartile range, i.e. it extends form the 25th to the 75th per- 
centile of data. The whiskers extend from the box until the ex- 
treme data. For each parameter the labels low, std and high 
refer to the values 1.00 - Ap t , 1.00, and 1.00 + Ap h 

The regression model coefficients, along with their statistical 
significance, are listed in Tables[3]and|4] In the first two columns 
of the Tables we report the least-squares estimates of the regres- 
sion coefficients and their errors; in the third column we report 
the f-statistic for the tests of the statistical significance of the 
covariates. In the last row of the Tables we report the residual 
standard error of the fit <x and the value of the squared multiple 
correlation coefficient R 2 , i.e. the fraction of the variance of the 
data explained by the model, defined as: 



R 2 _ j _ Wi-yi) 



(2) 



where y, are the values of the dependent variable, y their mean 
value, and yi the values predicted by the linear model. 

All the tests reach a very high statistical significance (p- 
values < 2 x 10~ 16 ). The physical relevance of the different co- 
variates can be assessed by looking at Fig. where we show 
the boxplots highlighting the influence of the variation of the 
inputs among the three values chosen for the calculations. For 
each parameter /?,• the boxplots are a convenient way to summa- 
rize the variability of the data subsetted according to the three 
values of p t (1.00 - A/?,-, 1.00, and 1.00 + A/?,, labeled as low, 
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std, and high in the plots). The black thick lines show the me- 
dian of the data set, while the box marks the interquartile range, 
i.e. it extends form the 25th to the 75th percentile of data. The 
whiskers extend from the box until the extreme data, the bottom 
whisker ranges from the sample minimum to the first quartile 
while the top whisker from the third quartile to the sample max- 
imum. While the position of the medians are related to the effect 
of the parameter in study in each plot, the extension of the box 
and whiskers are due to the variation of all the other parameters. 
The larger the separation of the medians with respect to the di- 
mension of the boxes and the greater the importance of a given 
parameter. 

The Figures show that the radiative opacity uncertainty 
largely dominates for both Lbto and ?h- For Lbto the impact 
of an increase of 5% in k r is AlogLBTo/Lo = 0.05 x (-0.276) = 
-0.0138 dex, while the impact of the variation of the second 
most important input - i.e. 14 N(p,y) I5 reaction rate - is only 
A log Lbto /Lq = -0.0028 dex. The sum of the effects in the sta- 
tistical models does not account for the whole variation of the 
data, due to the presence of the random variation extracted in the 
residuals of the statistical models. 

In the case of tu the dominance of k r is even larger since its 
variation accounts for A?h = 0.579 Gyr, and the second most 
important input, the microscopic diffusion velocities, for A?h = 
0.065 Gyr. 



Table 3. Fit of BTO log-luminosity (dex). 





Estimate 


Std. Error 


t value 


Impact 
(dex) 


Po 


6.00 x 10~' 


3.43 x 10~ 3 


174.91 




Pi (PP) 


6.76 x 10~ 2 


2.81 x 10~ 3 


24.06 


0.0020 


P2 ( 14 N) 


-2.80 x 10 2 


8.43 x 10"* 


-33.25 


-0.0028 


Pi (k) 


-2.76 x 10-' 


1.69 x 10~ 3 


-163.45 


-0.0138 


Pa (v d ) 


-9.56 x 10" 3 


5.62 x 10~ 4 


-17.00 


-0.0014 




cr = 6.2 x 10 


- 4 dex; R 1 = 


0.9974 





Notes. In the first two columns: least-squares estimates of the regression 
coefficients and their errors; third column: f-statistic for the tests of the 
statistical significance of the covariates. In the last column is reported 
the physical impact of the variation Ap, of the various inputs. All the 
p-values of the tests are < 2 x 10~ 16 . The residual standard error a and 
the squared multiple correlation coefficient R 2 are reported in the last 
row. 



In the cases of logL t i p , log Lhb and the linear models 
were: 

7 

log L tip , log L HB , M c He =fo + YjPiPi ( 3 ) 

where the three additional physical inputs considered are respec- 
tively the triple-o- reaction rate, the neutrino emission rate, and 
the conductive opacity k c (see Table [TJ. The regression model 
coefficients, along with their statistical significance, are listed in 
Tables [5] [6] and [7] The boxplots of the influence of the various 
inputs on the calculations are presented in Figs. [6] [7] and[8] 

For log L t ip and log Lhb the most important factor is again 
the radiative opacity k r , although its impact is not as dominant 
as in the previous cases. For instance, in the case of L t ; p , the 
impact of an increase of 5% in k r is AlogL t i p /LQ = -0.0149 
dex, while the variation of the second most important input, the 
triple-a reaction rate, accounts for A log L tiv /LQ = -0.0061 dex, 
i.e. about 40% of the k r effect. 

For Lhb the effect of the variation of k r is A log Lhb /Lq = 
-0.0216 dex, while the one of triple-c is A log Lhb /Lq = 
-0.01 15 dex, i.e. about 55% of the previous one. 

Regarding the helium core mass, the main variation is due to 
triple-or reaction rate: an increase of 20% of this value accounts 
for AM He = -0.00144 M Q ; the effect due to the uncertainty on 
neutrino emission rate, radiative and conductive opacities are all 
of the same order, i.e. about 45% of the triple-o- one each, while 
I4 N(p,y) ls O reaction rate uncertainty accounts for about 35% of 
the triple-a effect. 



Table 5. Fit of RGB tip log-luminosity (dex). 







Estimate 


Std. Error 


; value 


Impact 
(dex) 


Po 




3.68 


6.41 x 10~ 4 


5743.68 




Pi 


(PP) 


-2.56 x 10- 3 


4.11 x 10~ 4 


-6.23* 


0.0000 


Pi 


( 14 N) 


2.74 x 10- 2 


1.23 x 10~ 4 


222.28 


0.0027 


Pi 




-2.98 x 10~' 


2.47 x 10~ 4 


-1209.61 


-0.0149 


Pi 


(Vd) 


1.72 x 10~ 3 


8.22 x 10~ 5 


20.91 


0.0002 


Pi 


(3a) 


-3.05 x 10~ 2 


6.16 x 10~ 5 


-494.01 


-0.0061 


p 6 


(v) 


8.22 x 10~ 2 


3.08 x 10 4 


266.77 


0.0033 


Pi 


(k c ) 


-5.81 x 10~ 2 


2.47 x 10- 4 


-235.71 


-0.0029 






or = 4.7 x 10~ 4 dex; R 1 = 


0.9988 





Notes. The p-values of the tests are < 2 x 10 16 , if not differently 
specified. The column legend is the same as in Table[3] 
w p-value = 5.6 x lO" 10 . 



Table 4. Fit of central hydrogen exhaustion time (Gyr). 





Estimate 


Std. Error 


t value 


Impact 










(Gyr) 


Po 


-2.29 


0.0421 


-54.28 




Pi (PP) 


1.81 


0.0345 


52.46 


0.054 


P2 ( 14 N) 


-0.139 


0.0104 


-13.39 


-0.014 


Pi (kr) 


11.6 


0.0207 


559.81 


0.579 


Pa (v d ) 


-0.432 


0.00690 


-62.56 


-0.065 



a = 0.0076 Gyr; R 1 = 0.9998 



Notes. All the p-values of the tests are < 2 x 10 16 . The column legend 
is the same as in Table [3] 



The statistical models presented in this Section rely on the 
assumption of the linearity of the output of the calculations with 
respect to the perturbations of the physical inputs. As a check of 
this hypothesis, we made some additional runs varying only one 
parameter at a time, in a set of 9 points spanning the assumed 
range of uncertainty. The procedure is fully justified by the lack 
of interaction among physical inputs in the linear models pre- 
sented in Tables [3] - [7] In fact in presence of a relevant interac- 
tion, the values of the multiple correlation coefficient R 2 of the 
fits would be much smaller than the ones found in the models, 
which had a minimum value of 0.997 for Lbto, i-e- the 99.7% of 
the variance of the data are explained by the model. 

In Fig. |9]we present a graphic test of linearity. The four pan- 
els show, for each evolutionary feature studied in the Section, 
the results of the calculations for the variation of the two most 
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Table 6. Fit of helium core mass M" e (M©). 



Estimate 



Std. Error 



t value 



Impact 

(Mq) 



Po 




4.78 x 10- 


-i 


1.03 x 10- 


4 


4635.03 




Pi 


(PP) 


-1.20 x 10~ 


-3 


6.61 x 10- 


-5 


-18.16 


0.00004 


Pz 


( I4 N) 


-4.95 x 10 


-3 


1.98 x 10- 


-5 


-249.69 


-0.00049 


Ih 


(h) 


1.31 x 10" 


-2 


3.96 x 10- 


-5 


329.33 


0.00065 


Pa 


(Vd) 


1.93 x 10- 


-3 


1.32 x 10- 


-5 


146.35 


0.00029 


ft 


(3a) 


-7.20 x 10- 


-3 


9.91 x 10- 


-6 


-726.95 


-0.00144 


P 6 


(v) 


1.61 x 10- 


-2 


4.96 x 10- 


5 


324.28 


0.00064 


Pi 


(h) 


-1.20 x 10- 


-2 


3.96 x 10- 


5 


-303.07 


-0.00060 



o- = 7.6 x IP" 3 Mq- R 1 = 0.9976 



Notes. All the p-values of the tests are < 2 x 10 16 . The column legend 
is the same as in Table|3] 



Table 7. Fit of the ZAHB log-luminosity log L H b (dex) at log T e s 
= 3.83. 



Estimate 



Std. Error 



t value 



Impact 
(dex) 



Po 




2.03 


1.20 x 10- 


-3 


1696.88 




Pi 


(PP) 


6.93 x 10~ 3 


7.68 x 10- 


-4 


9.03 


0.0002 


P2 


( ,4 N) 


1.20 x 10~ 2 


2.30 x 10" 


-4 


52.04 


0.0012 


Pi 


(k) 


-4.33 x 10" 1 


4.61 x 10' 


-4 


-939.64 


-0.0216 


Pa 


(Vd) 


-1.12 x 10~ 2 


1.54 x 10- 


-4 


-73.17 


-0.0017 


Ps 


(3a) 


-5.78 x 10~ 2 


1.15 x 10- 


-4 


-501.83 


-0.0115 


P 6 


(v) 


5.77 x 10~ 2 


5.76 x 10" 


-4 


100.22 


0.0023 


Pi 


(h) 


-4.36 x 10~ 2 


4.61 x 10 


-4 


-94.70 


-0.0022 






a = 8.8 x 10 


dex; R 1 




0.9981 





Notes. All the p-values of the tests are < 2 x 10 16 . The column legend 
is the same as in Table|3] 



important physical inputs with superimposed the linear fit of the 
different trends. In all the cases the assumption of linearity is 
respected with high accuracy. 

Since the hypotheses of linearity in the inputs and their in- 
dependence hold, the statistical models presented in this Section 
can be safely used to interpolate the effect of a perturbation of 
the physical inputs in the assumed range. 

In the case of the central hydrogen exhaustion (i.e. X c = 
0), the effect of the radiative opacity variation will equate the 
one due to the second most important contribution, namely that 
due to the microscopic velocity variation, for a perturbation 
estimated as the product of assumed perturbation by the ratio 
of the effect of the two inputs as resulted by the model, i.e. 
0.05 x 0.065/0.58 = 0.0056. This means that in order to reduce 
the uncertainty in the predicted main sequence lifetime due to the 
radiative opacity at the level of the second largest contribution, 
the uncertainty affecting the Rosseland mean opacity should de- 
crease from the currently adopted 5% down to 0.56%. For the 
turn-off Lbto and ZAHB Lhb luminosities, the radiative opac- 
ity remains the most important uncertainty source for assumed 
perturbation greater or equal to 0.05 X 0.0028/0.014 = 0.01 and 
0.05 x 0.012/0.022 = 0.027, respectively. Thus, in order to re- 
duce the uncertainty in the turn-off luminosity caused by the ra- 
diative opacity to the same order of that caused by the second 
most important contribution, i.e. the 14 N(p,y) I5 reaction rate, 
the uncertainty of the former should be decreased from 5% to 
1%. Finally, the radiative opacity should be known with a preci- 
sion better than about 2.7% to lead a variation in the predicted 



-j 

CO 

o 

,g 

23 

o 



►J 

CO 

a 
,g 

23 

o 

OS 



low std high 
N (p, y) 15 



low 



std 
3a 



high 



-J 

CO 

o 

.Qi 
23 

a 



►J 
-J 

CO 

a 
.gl 

23 

o 

2< 



low std high 
Neutrinos 



low 



std 



high 



o 
,g 

23 

a 



-J 

50 

o 
,g 

23 

o 

2< 



low 



std 



high 



low std high 
Diffusion 



Fig. 6. Boxplot of the impact of the variation of the most impor- 
tant physical inputs on L t i p . 

ZAHB luminosity lower than the second largest contribution, i.e. 
the triple-a reaction rate. 

5. Global physical uncertainty in stellar isochrones 

The direct theoretical counterparts of the observed color- 
magnitude diagrams of simple stellar populations are isochrones 
rather than tracks with fixed mass. For this reason, we extended 
the previous statistical analysis of the physical uncertainties af- 
fecting stellar models to isochrones with age in the range 8-14 
Gyr, suitable for galactic globular clusters, with time steps of 1 .0 
Gyr. 

Since we were interested mainly in studying the variation 
near the turn-off region, we performed calculations varying 
only the four physical inputs p\, ... ,p\ that can influence the 
evolution until this phase, as shown in the previous section. 
Thus, for each set of perturbed input physics, i.e. for each set 
of multipliers values pi,... ,p4, we computed a grid of stellar 
tracks with different masses, for a total of 3 4 = 81 grids. Each 
grid contains 12 stellar models with mass values chosen to 
accurately reconstruct the zone near the BTO, namely M = 
0.40, 0.50, 0.60, 0.70, 0.80, 0.85, 0.87, 0.90, 0.92, 0.95, 1 .00, 1 . 10 
Mq. Then we computed 972 stellar tracks and 567 isochrones 
with fixed chemical composition (Z = 0.006, Y = 0.26) and 
mixing-length parameter (ff m i = 1.90). 

The global effect of the combined variation of the four 
selected inputs on the 12.0 Gyr isochrone in the (logr e ff, 
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0.90 0.95 1.00 1.05 1.10 0.85 0.90 0.95 1.00 1.05 1.10 1.15 

perturbation perturbation 




0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 

perturbation perturbation 

Fig. 9. Study of the linearity of the outputs of the calculations with respect to the variation of the physical inputs in their range of 
uncertainty. The symbols are the results of the calculations, while the lines are the linear fit of each set of values. Only the two most 
important inputs are shown in each panel. 



logL/L(T)) plane is displayed in Fig. [10] The stripe around the 
isochrone computed with standard inputs represents the region 
of the plane spanned by the perturbed stellar models. 

Figure QT| shows the boxplots of three selected quantities as 
a function of the isochrone age, namely the turn-off luminosity 
(Lg° Q ) and mass (Mg° ), and the differences between log Lhb 
and logLg° Q . The last one is the theoretical counterpart of the 
AV(TO-HB), i.e. the visual magnitude difference between turn- 
off and horizontal branch regions, used as age indicator in the 
vertical method technique. 

Table [8] lists the half-width of the variation range in predic- 
tions of the selected quantities for our reference stellar isochrone 
of 12 Gyr due to current input physics uncertainties. The turn-off 
log luminosity log Z^o and mass M^° Q vary of +0.013 dex and 
+0.015 Mq, respectively, while the log L HB /L l ^° vary of +0.05 
dex. 



Table 8. Range half-width of variation in theoretical predictions 
of selected quantities for our reference isochrone of 12 Gyr, with 
Z = 0.006 and Y = 0.26, due to input physics uncertainties. 



quantity 


range half-width 


lop / ,s " 
rug ^ BT0 


0.013 dex 


M is ° 
j "bto 


0.015 Mq 


log L„ s /L™ 


0.05 dex 



6. Statistical analysis of physical uncertainty in 
stellar isochrones 

To explore the influence of the physical inputs variation on the 
reconstructed turn-off log-luminosity log L^ Q and mass M l *° Q , 
we chose to pool the results for different ages removing the trend 
due to the age, increasing the possibility to detect the effect of 
the perturbations due to a larger statistic. The removal of the 
age trend was done by adapting a linear model to either log £g° 
and Mg° using as covariate a categorical variable, i.e. a variable 
which assumes different values for each age. The variable is then 
used as in a classical ANOVA (i.e. analysis of variance) model 
to filter out the variation due to the different isochrones ages. 
The residuals of the models were then regressed with respect to 
p\,...,Pi. The procedure guarantees the removal of the mean 
trend from the data, i.e. every set of residuals for each age has 
mean value of zero. The technique assumes that the effect of p, 
are the same at different ages in the studied time interval, i.e. 8- 
14 Gyr. As a rapid check of the hypothesis, we performed two 
Bartl ett tests of homogeneity of variances of residuals among 
ages (ISnedecor & Cochranl[l9 89). The tests did not suggest any 
problem (log L™ Q p-value = 0. 16, M*° Q />value = 0.26). 

The results of the linear models for log Z^o and M l ^° Q are 
in Table l9l and [TUl respectively. In Figs. [T2l and [T3l we present 
the boxplots which evidence the influence of the variation of the 
physical inputs among the three values chosen for the calcula- 
tions. For the turn-off mass M l *° Q , the effect of the variation of 
the radiative opacity is largely dominant. An increase of 5% in k r 
produces a variation of AM = 0.05 x 0.270 = 0.0135 Mq, while 
the impact of the variation of the second most important input 
- i.e. microscopic diffusion velocities - is only AM = -0.0014 
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Fig. 11. Left panel: boxplots for BTO luminosity due to the variation of the selected physical inputs as a function of the isochrone 
age in Gyr. Central panel: same plot for mass at BTO. Right panel: same plot for difference of HB and BTO log luminosity. 
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Fig. 7. Boxplot of the impact of the variation of the most impor- 
tant physical inputs on M^ e . 
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Fig. 8. Boxplot of the impact of the variation of the most impor- 
tant physical inputs on Lhb- 



Mq. The total range of variation of AM™ is [-0.0184 - 0.0180] 
M Q . 

For the luminosity of the turn-off L"° Q , no physical input 
definitely dominates on the others. The most important factor 
turned out to be radiative opacity with a variation A log L/Lq = 
0.0042 dex, followed by 'H(p,ve + ) 2 H reaction rate with a vari- 
ation A log L/Lq = 0.0031 dex, 14 N(p,y) 15 reaction rate with 
A log L/Lq = -0.0029 dex, and microscopic diffusion velocities 
with A log L/Lq = -0.0026 dex. The total range of variation of 
AlogL^ is [-0.0142 - 0.0135] dex, about 2/3 of the one found 



in Sec. [3] for the turn-off luminosity logLBTO of the reference 
track. This lower uncertainty in the isochrones with respect to 
the tracks used to build them is due to the fact that isochrones 
of the same age but computed with different sets of multipliers 
Pi,...,P4 values have turn-off regions populated by different 
masses M™° Q . As already shown, the variation of k r splits the 
mass at BTO in three separate sets for each value of the param- 
eter (see Fig. fT3l . The higher the value of k r , the more massive 
are the stellar models at BTO. The luminosity of these models, 
which are intrinsically higher, are depleted by the high value of 
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I Table 9. Fit of the isochrones BTO log-luminosity (dex). 



-j 



3.76 3.75 
log T ct , 



Fig. 10. HR diagram showing the error stripe due to the variation 
of the physical inputs for a 12.0 Gyr isochrone (zoom of the TO 
region). 



the radiative opacity, as in the third panel of first row in Fig. 
The opposite behavior was evidenced for low value of k r . The 
net effect is a shrinkage of the luminosity range spanned by the 
models. 

The most important and common application of theoretical 
isochrones is the age estimate of stellar clusters. As such, it 
would be of primary importance to evaluate the current uncer- 
tainty affecting the age inferred from the comparison between 
observed and predicted turn-off luminosity, in addition to the un- 
certainty in the BTO luminosity at a given age analyzed above. 
The computed isochrones can be used for such purpose, even 
if a special care is needed. In fact a direct modeling of the age 
with respect to log £gj is not feasible, since the age is used in 
the calculation as a parameter with null uncertainty and it can- 
not be viewed as random variable, as required for a linear model 
dependent variable. Therefore we proceed in the following way. 
We found that the relation between log L l ^° Q and the age is well 
described by the following linear model: 



logL; 



BTO 



(4) 



where tg = log t, with t the age of the isochrone in Gyr. The 
parameters of the model are: fio = 1.2871, fii = -0.8771, with 
cr = 0.0056 dex. Since we do include explicitly the contribu- 
tion of the parameter pi,.. -,pa, they contribute to the inflation 
of the residual standard error of the model, reflecting our global 
uncertainty. This model can be used to perform a reverse infer- 
ence on the value of the age, given L™° Q . As described in detail 
in Appendix|B] in the explored range [8 - 14] Gyr, a typical un- 
certainty in the age is about + 375 Myr. 

With an analysis like the one described in the previous para- 
graph for L™° Q , we can obtain a reverse inference of the uncer- 
tainty on age given the value of log Lhb / ^bto- ^h e linear model 
adapted to data is: 



log L HB IU 



BTO 



■■fa +/3i t 9 



(5) 



where tg = log t, with t the age of the isochrone in Gyr. The 
parameters of the model are: Pq = 0.2779, fi\ = 0.8771, with 
cr = 0.024 dex. As a result of the large value of cr, a typical 
uncertainty in age, given log Lhb/Lbto' * s °^ tne OT< ^ eT °f ± 1-25 
Gyr. 





Estimate 


Std. Error 


t value 


Impact 
(dex) 


Po 


-1.40 x 10-' 


2.72 x 10- J 


-51.51 




Pi (pp) 


1.04 x 10-' 


2.23 x 10- 3 


46.60 


0.0031 


Pi C 4 N) 


-2.95 x 10~ 2 


6.68 x 10- 4 


-44.13 


-0.0029 


Ps (kr) 


8.34 x 10~ 2 


1.34 x lO- 3 


62.43 


0.0042 


Pa (v d ) 


-1.77 x 10-2 


4.45 x 10~ 4 


-39.65 


-0.0026 


tr = 0.0013 dex; R 1 = 0.945 



Notes. All the p- values of the tests are < 2 x 10 16 . The column legend 
is the same as in Table[3] 



Table 10. Fit of mass at BTO (M ). 







Estimate 


Std. Error 


t value 


Impact 

(M Q ) 


Po 




-3.00 x 10- 1 


2.08 x W' 3 


-143.97 




Pi 


(PP) 


4.13 x 10~ 2 


1.71 x lO" 3 


24.23 


0.0012 


Pi 


( 14 N) 


-2.40 x 10~ 4 


5.12 x 10~ 4 


-4.69* 


-0.0002 


Pi 


(kr) 


2.70 x 10-' 


1.02 x 10 3 


264.03 


0.0135 


Pa 


(Vd) 


-9.31 x 10" 3 


3.41 x 10~ 4 


-27.30 


-0.0014 






cr = 9.9 x 1C 


Mq\ r 1 = 


0.992 





Notes. The p-values of the tests are < 2 x 10 16 , if not differently 
specified. The column legend is the same as in TablefJ] 
w p-value = 3.36 x 10~ 6 . 
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Fig. 12. Boxplot of the impact of the variation of the chosen 

3 turn-off luminosity L^° . The 
14] Gyr are pooled together, see 



physical inputs on the isochrone turn-off luminosity L™° Q . The 
results for ages in the range [8 
text for details. 



7. Conclusions 

In this paper we addressed the problem of a quantitative 
and systematic evaluation of the cumulative propagation 
of physical uncertainties in current generation of stellar 
models of low mass stars from the main sequence to the 
zero age horizontal bra nch. At varian c e with several pre- 
vious work (see e.g. IChabover et al.l Il995t ICassisi et al.l 
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Fig. 13. Boxplot of the impact of the variation of the chosen 
physical inputs on the isochrone turn-off mass Mg° Q . The re- 
sults for ages in the range [8 - 14] Gyr are pooled together, see 
text for details. 



1998 1 ; ICastellani & DegFInnocentil 119991; ICastellani et al l 



2000; llmbriani et all 1200 U iPrada Moroni & Stranierol 12002; 
Safaris et all l2002t llmbriani et all 12004 I Weiss et all I2OO5I: 
Prada Moroni & Stranierol2007l:IValle et al.l2009l:lTognelli et all 
201 ll) . where a single input physics was changed at a time, we 
performed a systematic and simultaneous variation on a fixed 
grid within their current range of uncertainty, in a way to obtain 
a full crossing of the perturbed input values, of the main physical 
inputs adopted in stellar codes (i.e. 'H(p,ve + ) 2 H, 14 N(p,y) 15 0, 
and triple-a reaction rates, radiative and conductive opacities, 
neutrino energy losses, and microscopic diffusion velocities). 

Although very expensive from the computational point of 
view, such an approach has the important advantage with respect 
to the previous one to be more robust against possible interac- 
tions among the varied input physics, as any a priori indepen- 
dence among them is assumed. 

Relying on a set of stellar models fully covering all the possi- 
ble combinations of simultaneously perturbed input physics, we 
were able to compute the error stripe associated to our reference 
stellar track, i.e. M = 0.90 Mq with initial metallicity Z = 0.006 
and helium abundance Y = 0.26, from the pre-main sequence 
up to the RGB tip, in different planes (i.e. logL vs log T e g, X c 
vs age, A logL vs model, and A log/? vs model). We built also 
the error stripe associated to the ZAHB locus resulting from the 
evolution of a progenitor of M — 0.90 Mq. As far as we know, 
this is the first time that an error stripe is computed and plotted 
for stellar tracks (see Appendix lAl for a detailed description of 
the error stripe construction method). 

Furthermore, we quantified the extension of the global vari- 
ability regions (i.e. due to the cumulative effect of all physi- 
cal uncertainties) for some relevant stellar quantities, without 
any assumption on the parent distributions of the varied phys- 
ical inputs. In particular, we focused on the turn-off luminos- 
ity Lbto, the central hydrogen exhaustion time tu, the luminos- 
ity L t ip and the helium core mass M^ e at the RGB tip, and the 
ZAHB luminosity in the RR Lyrae region Lhb, were computed, 



too (see e.g. Table|2|. We found that the turn-off log luminosity 
logLBTO varies of +0.021 dex, while the RGB tip logL t ; p and 
ZAHB log Lhb ones of +0.03 dex and +0.045 dex, respectively. 
Thus, uncertainties of the order of ±0.075 mag and ±0.10 mag 
should be taken into account in distance modulo estimates ob- 
tained by theoretically calibrated RGB tip and ZAHB luminosi- 
ties. The predicted central hydrogen exhaustion time ?h varies of 
±0.72 Gyr, whereas the helium core mass at the RGB tip M^ e of 
±0.0042 Mq. 

This large computational effort (i.e. 3 159 stellar tracks) made 
possible a thorough and rigorous statistical analysis of the effects 
of the variations of the quoted physical inputs on the chosen rel- 
evant stellar evolutionary features, allowing to disentangle the 
contributions of the different inputs physics. 

The results of our extended statistical analysis show that the 
radiative opacity is, by far, the dominant source of physical un- 
certainty in almost all of the examined stellar evolutionary fea- 
tures, with the exception of the helium core mass M^ e at the 
RGB tip which is affected mainly by the uncertainty in the triple- 
a reaction rate. As an example, for the turn-off log luminosity 
log Lbto, in order to reduce the impact of the radiative opacity 
variation to the level of the second most important input, i.e. the 
14 N(p,y) 15 reaction rate, the uncertainty of the former should 
decrease from the current 5% down to 1 %. For the central hydro- 
gen exhaustion time tu, the prevalence of the radiative opacity is 
even larger, as for reducing its effect at the level of the second 
most important input physics, i. e. the microscopic diffusion ve- 
locities, its uncertainty should be decreased at the 0.56% level. 
Notice that the 5% uncertainty in radiative opacity assumed in 
our calculations is probably an underestimate, at least in some 
regions of the temperature-density plane of interest for stellar 
models. 

Beside this analysis of the cumulative uncertainty in a stel- 
lar track of fixed mass and chemical composition, we extended 
the previous analysis of the cumulative uncertainty due to the 
main input physics to theoretical isochrones in the age range 8- 
14 Gyr. We focused in particular on the turn off luminosity Lg° 
and mass M'^° Q , and log Lhb / log Lg° , which is the theoreti- 
cal counterpart of the visual magnitude difference between turn- 
off and horizontal branch regions (i.e. AV(TO - HE) used as 
age indicator in the "vertical method"). For an age of 12 Gyr, 
we found that the isochrone turn-off log luminosity log Lg° 
varies of ±0.013 dex, about 2/3 of the value found for the turn- 
off luminosity Lbto of the reference stellar track. For the same 
age, the mass at the isochrone turn-off varies of ±0.015 Mq 
and log Lhb / log L™ Q of ±0.05 dex. The large set of perturbed 
isochrones allowed us to perform the same kind of statistical 
analysis previously applied to tracks of fixed mass, with the aim 
to evaluate the effects of the single varied input physics. 

Finally, the availability of isochrones computed by varying 
simultaneously all the main input physics made possible to eval- 
uate the physical uncertainty affecting the age inferred from two 
of the most used cosmic clocks, namely the turn-off luminosity 
and the vertical method. For given Lg° and log Lhb / log Lg^ 
values, the inferred age varies in a of about ± 0.375 Gyr and ± 
1 .25 Gyr, respectively. 

An equally systematic approach is provided by Monte 
Carlo si mulations as those presented by Chaboyer and collab- 
orators dChabover et al.lll996l Il998b IChabover & Kraussl 120021: 



iKrauss & Chaboyer; 120031: iBiork & Chaboverll2006l) . However, 
a direct comparison with their results is difficult since the stud- 
ies address different questions about the global uncertainty. The 
main problem is due to the fact that we focus on the uncertainty 
of evolutionary features different from the one studied in those 
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papers. The only possible feature we can compare is the change 
in the evolutionary age due to the variation of radiative opacity. 
In IChabover & Kraussl d2002l) a 2% change in radiative opacity 
accounts for a 2.6% variation (i.e. 6.5% for a variation of ra- 
diative opacity of 5%) of the age of a Sub-Giant Branch star. 
As proxy of this quantity, we found that a 5% increase in radia- 
tive opacity account for a change of the time of central hydro- 
gen exhaustion of 5.5%, i.e. 1 % lower to the equivalent variation 
quoted above. The small difference can be in principle due to the 
diffe rent sampling schema for the radiative opacity adopted in 
Chabover & Krauss (120021) . and to the different chemical inputs 
and solar mixture employed in the calculations. 

Although the present paper addresses some fundamental top- 
ics about the uncertainty of modern stellar models, many other 
questions remain open and need further investigations. The first 
problem concerns the possibility to extend the results of our 
analysis - performed at fixed values of metallicity and initial 
helium abundance - to different values of the chemical inputs. 
In the lack of the very huge set of computations required to 
specifically address this topic, the extrapolation of the results 
presented in this paper to values of Z and Y very different from 
the ones employed in the computations is to be considered with 
care. Some cautions should also be adopted for "borderline" stel- 
lar models which can develop a different structure - such as the 
development of a convective core - due to a perturbation of a 
physical input. In these few cases the evaluation of the uncer- 
tainty can not be inferred but has to be performed with direct 
computations. 

A second question concerns the possibility to extend the re- 
sults presented here to the computations performed with other 
stellar evolutionary codes. This point could be addressed only 
with the replication with other codes, with the very same in- 
puts and configurations, of the computations presented here. The 
availability of such calculations by different groups of the stellar 
evolution community would be of invaluable importance. In fact 
in this way could be assessed not only the uncertainty due to the 
variations of the physical inputs, but also to quantify the impact 
of another source of uncertainty, i.e. the possible systematic dif- 
ference among various codes due to different implementations 
of physical mechanisms and to the algorithms used in the com- 
putations. 
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Appendix A: Tracks reduction and error stripe 
construction on the theoretical plane 

The first step of the construction of the stellar track uncertainty 
stripe on the theoretical plane is the reduction of the raw tracks 
to a set of tracks with the same number n of homologous points. 
The reduction is based upon the identification of some keystone 
evolutionary points on the raw track and subsequent deployment 
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of the same number of points - obtained by interpolation - on 
all the tracks. 

As reference points we adopted the following: 

1 . PMS 1 : the point for which gravitational luminosity reaches 
0.996 times the surface luminosity. 

2. PMS2: the point for which the gravitational luminosity de- 
creases of 0.05 from the value at PMS1. 

3. ZAMS: the point for which the central hydrogen abundance 
drops below 99% of its initial value. 

4. MSI: the point for which the central hydrogen abundance 
drops below 7%. 

5. MS2: the point for which the central hydrogen abundance 
drops below 1%. 

6. MS 3: the point for which the central hydrogen abundance 
drops below 0.1%. 

7. HC: the point for which the central hydrogen is exhausted. 

8. RGB1, or RGB start: the point on the track at maximum dis- 
tance from the line connecting the HC and the RGB2 point. 

9. RGB2, or RGB bump: the point for which the luminosity of 
RGB start decreasing. 

10. RGB3, or He-flash: the point for which the He burning lumi- 
nosity reaches 100 times the surface luminosity. 

The ZAHB models have been computed through a synthetic 
method where the starting models were obtained by accreting 
envelopes of different mass extensions onto the He-core left at 
the tip of the RGB. Then, full evolutionary calculations were 
started again as thermal relaxed models in the central He burn- 
ing phase; ZAHB point corresponds to the model in which the 
equilibrium abundance of CNO burning secondary elements is 
reached, after about 1 Myr. 

Let {T,}, i = 1, . . . , N be the set of the N reduced tracks and 
let Ti(j) be the j'th point on the ith reduced track. Let's define: 



T{j) = {T i {j)} i = l, 



(A.l) 



the set of the jth points over the whole set of reduced tracks. Let 
H C {X) be the convex hull of a generic set X. We define: 



(k+\ 



<H(k) = H c 



{JT(J) 

V j=k 



(A.2) 



the convex hull of the set composed by the A:th and (k + l)th 
points of the reduced tracks. 

The full stripe is then given by: 



(A.3) 



k=\ 



The required computations we re carried out using R 2.14.1 
dR Dev elopmen t Core Teamll201 ll) . 



Appendix B: Reverse inference on the age of a 
isochrone given the TO luminosity 

Let us consider a simple linear model: 

y =v8 +A x (B.l) 

the 95% confidenc e interval (yi, V 2) on a future observation y 
given x, is (see e.g. Faraway 2004): 



i_„/9 I 1 (x — x) 2 
yia =p 0+ p l x±tlfcrJl + - + { -^ 



where Bq, B\, & are the least-squares estimates of the model pa- 
rameters, a is the required confidence level (in our case, a = 
0.05), n is the number of points in the model, t 1 ^ 2 is the 1 -a/2 
quantile of the Student t distribution with n - 2 degrees of free- 
dom, x is the sample mean value of x, d 2 is the sample deviance 
of x. 

In our case we have: y = logLg° and x = logT, with 
T the age of the isochrone in Gyr. The boundaries of the 95% 
confidence interval are displayed along with the best fit line in 
Fig. IB. II where we show the construction of the range of un- 
certainty on age, given the value of BTO log luminosity. As an 
example, for log L™° Q = 0.35 dex the estimated ages from the 
models lie in the range [1 1.37 - 12.05] Gyr. 

Appendix C: On the mixing-length influence 

The results quoted in the present paper are computed for a fixed 
value of mixing-length parameter, i.e. a m i = 1.90. The possibil- 
ity to extend them to different values of this parameter must be 
checked by computing stellar models with different a m \. 

In presence of an effect due to the mixing-length, we ex- 
pect a difference in the regression coefficients for the models 
summarized in Tables |3]|7] when they are computed from stel- 
lar model with different a m \. The computation of the huge sets 
of stellar models for different mixing-length values is not fully 
needed since in the linear models presented in Sec. [4] is shown 
that the various physical inputs do not interact. Therefore a sub- 
set of stellar models carefully selected can suffice to assess the 
presence of a mixing-length effect of distortion on the regression 
coefficients. 

To perform the analysis we computed, for the two values of 
the mixing-length parameter a m \ = 1.70, 1.80 (a m i = 1.74 is the 
solar-calibrated value of the mixing-length parameter), 8 1 stel- 
lar models each. The models to be computed were randomly 
selected using a latin hypercube sampling design, which is an 
extension of latin square to higher dimensions and has optimal 
property in reducin g the varian ce of the estimators obtained from 
the linear models dSteinl 1 1 9871) . The random selection was per- 
formed using the R library lhs dCarneilll2012l) . The correspond- 
ing 81 models for a m \ - 1.90 were extracted from the 2187 orig- 
inal calculations and these 243 stellar models are used for estab- 
lishing the influence of a m i. 

The analysis is performed by adapting to data linear models 
of the form: 



logL B TO, ?H = A) + i^jfii Pi) * a ml 



i=l 



or: 



logL t ip, logL HB , =/3 + C^Pi pd * a mX 



(CI) 



(C2) 



(B.2) 



where, for convenience, we adopted the operator "*" defined as 
A*B=A + B+ AB. The effect of distortion of the regression 
coefficients due to the presence of the mixing-length parameter 
is represented by the interaction between each multiplier pi and 
a m \. The statistical significance of these coefficients is of lim- 
ited importance, since it can be arbitrarily increased by select- 
ing a larger subsample thus reducing the standard error estimate 
of the coefficients. A more useful indicator is the physical im- 
pact of the introduction of the interaction in the linear models. 
As explained in the text the impact of a perturbation Ap, on a 
parameter p-, which enters a linear model can be evaluated as 



15 



Valle, G. et al.: Cumulative physical uncertainty in modern stellar models 




8 9 10 11 12 13 14 8 9 10 11 12 13 14 



age (Gyr) age (Gyr) 

Fig. B.l. Left panel: reverse inference on the isochrone age given the BTO log luminosity. The solid line represents the best fit 
model, the dashed lines define the 95% confidence interval on the prediction of the model. The dotted lines show the uncertainty on 
age given a value of BTO log luminosity. Right panel: same for log LuB/Lgj . 



Api ■ Pi, where /?, is the estimate of the regression coefficient. 
In presence of the interaction the impact of the same perturba- 
tion Api combined with a variation of Aa m \ can be estimated as 
Apt ■ Aa m \ ■ fij where fij is the regression coefficient of the inter- 
action term under analysis. In table IC.ll we report the impact of 
the interaction terms in the various models for Aa m \ = 0. 1 . In all 
cases the physical impact of the interaction is negligible, so we 
can conclude that the regression presented in Sec. [4] are robust 
for an acceptable change in a m \. 



Table C.l. Additional impact of the variation of the mixing- 
length value by 0.1 on the physical impacts presented in Sec. [4] 



source 


log L BT0 


to 


log L tip 


M L Hc 


logL HB 




(dex) 


(Gyr) 


(dex) 


(M©) 


(dex) 


»ml • PP 


0.00002 


0.0008 


-0.00006 


1.3 x 10- & 


0.00004 


aw 14 N 


0.00002 


0.0000 


0.00002 


2.1 x 10- 7 


-0.00001 


«ml • K 


0.00026 


0.0033 


0.00002 


-6.0 x 10~ 6 


-0.00006 


«ml • V d 


0.00000 


0.0001 


0.00000 


2.7 x 10~ 7 


0.00001 


a m] ■ 3a 






-0.00016 


-4.0 x 10~ 6 


-0.00009 


«ml • V 






0.00014 


1.7 x 10~ 6 


0.00002 


Oral • k c 






-0.00006 


-5.4 x 10~ 7 


-0.00008 



Table D.l. Impact of equation of state change from OPAL to 
FreeEOS for the examined evolutionary features. 



Evolutionary feature 


OPAL 


FreeEOS 


Difference 


logL BTO (dex) 


0.3548 


0.3520 


0.0028 


f H (Gyr) 


10.557 


10.461 


0.096 


log L tlp (dex) 


3.4055 


3.4075 


-0.0020 


M L Hc (Mq) 


0.4835 


0.4850 


-0.0015 


log L HB (dex) 


1.5634 


1.5704 


-0.0070 



Notes. First column: evolutionary stage; second column: values ob- 
tained with OPAL EOS; third column: values obtained with FreeEOS; 
last column: differences among the two values. 

the impact of the EOS change. By comparing them with the val- 
ues in the last column of Table |2 one sees that the EOS change 
accounts for about 1/7 - 1/8 of the total variation due to all the 
other physical inputs. The only exception is the He core mass, 
for which the EOS change accounts for 1/3 of the total variation. 

Unlike the other analysis performed in the paper, in this case 
it is not possible to identify which of the varied thermodynal 
quantities mainly affects the total result. One also should be 
aware that a given thermodynamical quantity can have different 
influences in different evolutionary phases. 



Appendix D: On the EOS influence 

The assessment of the EOS influence on the stellar models can 
not be done in the same way of the other physical inputs, given 
the available information. As discussed in the text, the thermody- 
namic quantities required for computing stellar models and pro- 
vided by EOS are related together in a not trivial way, hamper- 
ing a simple parametrization of the uncertainty associated with 
them. 

Nevertheless, a rough estimate of the impact due to the 
present uncertainty on EOS can be computed in an alternative 
way. For this purpose, in this Appendix we provide the output 
of the stellar models computed with two different and widely 
adopt ed cho ices for equation-of-state: OPAL EOS and FreeEOS 
dlrwinll2004l ). 

Table IDTI reports - for each evolutionary feature - the val- 
ues obtained using the two different EOS and, in the last column, 
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